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Abstract. Faraday rotation of the polarization plane in magnetized thermal plasma provides one of the most 
efficient methods to deduce regular magnetic fields from radio astronomical observations. Since the Faraday 
rotation measure RM is proportional to an integral, along the line of sight, of magnetic field weighted with 
thermal electron density, RM is believed to yield the regular magnetic field averaged over large volume. Here we 
show that this is not the case in a turbulent medium where fluctuations in magnetic field and electron density are 
not statistically independent, and so contribute to RM. For example, in the case of pressure equilibrium, magnetic 
field can be anticorrelated with plasma density to produce a negative contribution. As a result, the strength of 
the regular magnetic field obtained from RM can be underestimated if the fiuctuations in electron density and 
magnetic field are neglected. The anticorrelation also reduces the standard deviation of RM. We further discuss 
the effect of the positive correlations where the standard treatment of RM leads to an overestimated magnetic 
field. Because of the anisotropy of the turbulent magnetic field, the regular magnetic fields strength, obtained from 
synchrotron emission using standard formulae, can be overestimated. A positive correlation between cosmic-ray 
number density and magnetic field leads to an overestimate of the strengths of the regular and total fields. These 
effects can explain the difference between the strengths of the regular Galactic magnetic field as indicated by 
RM and synchrotron emissivity data and reconcile the magnetic field strength in the Solar vicinity with typical 
strength of regular magnetic fields in external galaxies. 
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1. Introduction 

Estimates of magnetic field strength in the diffuse inter- 
stellar medium (ISM) of the Milky Way and other galax- 
ies are most efficiently obtained from the intensity and 
Faraday rotation of synchrotron emission. Other methods 
are only sensitive to relatively strong magnetic fields that 
occur in dense clouds (Zeeman splitting) or are difhcult to 
quantify (optical polarization of star light by dust grains) . 
The total / and polarized P synchrotron intensities and 
the Faraday rotation measure RM are integrals over the 
path length L, so they provide a measure of the average 
magnetic field in the emitting or magneto-active volume: 



where n^r and ric are number densities of relativistic and 
thermal electrons, B is the total magnetic field compris- 
ing a regular B and random b parts, B = B + b with 
(B) = B, (6) = and (B^) = # + {b'^), angular brackets 
denote averaging, subscripts _L and || refer to magnetic 
field components perpendicular and parallel to the line of 
sight, and K and Ki = 0.81 radm^^cm'^ /iG^^ pc^^ are 
certain dimensional constants. The degree of polarization 
p is related to the ratio (b'^) / B , 



P 

P = 7 



(2) 



I = K I n„Bl ds , 



P = K I ncrB]_ ds , 



RM = ^1 J ^cB\i ds , 

Send offprint requests to: R. Beck 



where the random field b has been assumed to be isotropic, 
TT-cr is assumed to be a constant, and po ~ 0.75 weakly de- 
pends on the spectral index of the emission (Burn 119661 
(1) Sokoloff et al. I1998|l . This is an approximate relation. In 
particular, it does not allow for any anisotropy of the ran- 
dom magnetic field (see Sect. 14. 1|) . for equipartition be- 
tween magnetic fields and cosmic rays (see Sect. I4.2|l and 
for depolarization effects; some generalizations are dis- 
cussed by SokoloflF et al. (|1998|). 
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Since ricr is difficult to measure, it is most often as- 
sumed tliat magnetic fields and cosmic rays are in energy 
equipartition; this allows one to express ricr in terms of B. 
The physical basis of this assumption is the fact that cos- 
mic rays are confined by magnetic fields. An additional as- 
sumption involved is that the energy density of relativistic 
electrons responsible for synchrotron emission (and so ricr) 
is proportional to the total energy density of cosmic rays; 
it is usually assumed that the energy density of relativistic 
electrons in the relevant energy range (i.e. several GeV) 
is one percent of the proton energy density in the same 
energy interval (Chapter 19 in Longair ll994|l . The valid- 
ity of this assumption requires that the diffusion/escape 
losses of the cosmic-ray electrons dominate over radiative 
losses, so that the energy spectrum of the electrons is not 
steeper than that of the protons fBeck [TM7|l . 

Estimates of B in our Galaxy via ricr determined from 
7-ray emission (Strong et al. I2000|l are generally consis- 
tent with the equipartition values (E. M. Berkhuijsen, in 
Beck|2ninj- However, Eq. Q is not consistent with the 
equipartition or pressure balance between cosmic rays and 
magnetic fields insofar as it assumes that rid- — const. 
Therefore, the regular magnetic field strength obtained 
using Eq. ^ can be inaccurate (see Section 0]). 

The thermal electron density ric in the ISM can be ob- 
tained from emission measure EM oc ds, although 
this involves additional assumptions regarding the filling 
factor of interstellar clouds. In the Milky Way, the dis- 
persion measures of pulsars, DM = ric ds provide in- 
formation about the mean thermal electron density, but 
the accuracy is limited by our uncertain knowledge of dis- 
tances to pulsars. Estimates of the strength of the regular 
magnetic field in the Milky Way are often obtained from 
the Faraday rotation measures of pulsars simply as 



This estimate is meaningful if magnetic field and thermal 
electron density are uncorrelated. If the fiuctuations in 
magnetic field and thermal electron density are correlated 
with each other, they will contribute positively to RM and 
Eq. Q will yield overestimated In the case of anti- 
correlated fluctuations, their contribution is negative and 
Eq. ||2Jl is an underestimate. In order to quantify this ef- 
fect, one needs a suitable model for the relation between 
magnetic fields and thermal electron density. As we show 
in Section 121 physically reasonable assumptions about the 
statistical relation between magnetic field strength and 
electron density can lead to Eq. Q being in error by a 
factor of 2-3 even in a statistically homogeneous magneto- 
ionic medium. Lerche fl9705 discussed the effects of cor- 
related fluctuations in magnetic field and electron density 
on Faraday rotation measures. Some results of that paper 
are presented in a questionable form, although the general 
conclusion agrees with that proposed here. 

Equation (PJ can also lead to significantly underesti- 
mated strength of the regular magnetic field if it is en- 
hanced in the interarm regions whereas electron density is 



maximum within the arms, as in galaxies with magnetic 
arms (BeckEOHIll- 

2. Magnetic field estimates 

The observable quantities |^ provide extensive data on 
magnetic field strengths in both the Milky Way and ex- 
ternal galaxies (Ruzmaikin et al. 119881 Beck et al. 119961 
Beck TMM IM1T|) . The average total field strengths in 
nearby spiral galaxies obtained from total synchrotron 
intensity / range from B ~ 4/iG in the galaxy M31 
to ~ 15 /iG in M51, with the mean for the sample of 
74 galaxies of B ~ 9 ± 3/xG (Beck 2Q0QJ. The typi- 
cal degree of polarization of synchrotron emission from 
galaxies at short radio wavelengths is p = 10%-20%, so 
Eq. Q gives B/B — 0.4-0.5; these are always lower limits 
due to limited resolution of the observations. Most exist- 
ing polarization surveys of synchrotron emission from the 
Milky Way, having much better spatial resolution, suf- 
fer from Faraday depolarization effects and missing large- 
scale emission and cannot provide reliable values for p (see 
Sect. EJ. Phillipps et al. (1981) obtained B/B = 0.6-0.7 
from analysis of the total synchrotron emission from the 
Milky Way along and perpendicular to the spiral arms. 
Heiles II1996|I derived similar values from starlight po- 
larization data. The total equipartition magnetic field in 
the Solar neighbourhood is estimated as = 6 ± 2 fiG 
from the synchrotron intensity of the large-scale, diffuse 
Galactic radio background (E. M. Berkhuijsen, in Beck 
IM7T|) . Combined with B/B = 0.65, this yields a strength 
of the local regular field of i? = 4± 1 fiG. Hence, the typi- 
cal strength of the local Galactic random magnetic fields, 
(fe2)i/2 ^ (^2 _ b'^)1/2 ^ 5 _^ 2_^G, exceeds that of the 
regular field by a factor {b'^y/'^/B = 1.3 ± 0.6. RM data 
yield similar values for this ratio (Sect. IV. 4 in Ruzmaikin 
et al.CnEHI; Ohno & Shibata[^p. 

Meanwhile, the values of B in the Milky Way ob- 
tained from Faraday rotation measures seem to be sys- 
tematically lower than the above values. RM of pulsars 
and extragalactic radio sources yield B = 1.4 ± 0.3 /xG 
in the local (Orion) (Rand & Lyne 119941 Frick et al. 
I2001|l . B = 1.7 ± 0.3 /iG in the Sagittarius-Carina spi- 
ral arm (Frick et al. 2001) ; for the Perseus arm, Frick et 
al. (EDOl) obtained B = 1.4 ± 1.2 ^G, and Mitra et al. 
(UnnSI), ^ = 1.7± 1.0 nG. The median value of S is about 
twice smaller than that inferred from other methods. 

There can be several reasons for the discrepancy be- 
tween the estimates of the regular magnetic field strength 
from Faraday rotation and synchrotron intensity. Both 
methods suffer from systematic errors due to our uncer- 
tain knowledge of thermal and relativistic electron den- 
sities, so one cannot be sure if the difference is signif- 
icant. Nevertheless, the discrepancy seems to be worry- 
ing enough to consider carefully its possible reasons. (We 
should emphasize that the main results of this paper are 
independent of whether or not the discrepancy is real.) 

The discrepancy can be explained, at least in part, if 
the methods described above sample different volumes. 



R. Beck et al.: Systematic bias in interstellar magnetic field estimates 



3 



The observation depth of total synchrotron emission, 
starlight polarization and of Faraday rotation measures 
are all of the order of a few kpc. Polarized emission, how- 
ever, may emerge from more nearby regions (see Sect.|Sl). 
However, there are more fundamental reasons for the dis- 
crepancy that we discuss in what follows. 

3. The effects of magneto-ionic fluctuations on 
RM 

In this section we show that a statistical correlation be- 
tween electron density fluctuations and turbulent mag- 
netic fields can affect significantly regular magnetic field 
estimates obtained from Faraday rotation measures. The 
effect of fiuctuations in the magneto-ionic medium on 
RM can be the main reason for the discrepancy between 
magnetic field estimates discussed in Sect. |3 Note that 
the effect discussed here arises from correlations at scales 
^ 100 pc (i.e., less than the basic turbulent scale), rather 
than from any relations between the averaged quantities 
(e.g., arising from the vertical stratification of the ISM or 
galactic density waves). 

In order to quantify the effect, one needs a specific 
physical model for the connection between thermal gas 
density and magnetic field. An appealing idea is based on 
the assumption that the ISM is in pressure balance involv- 
ing not only thermal pressure (Field et al. 119691 McKee 
& Ostriker [l977|l but also turbulent, magnetic and cosmic 
ray pressures (e.g., Parker IT979I Boulares & Cox 119901 
Fletcher & Shukurov l2()()l|l . Pressure balance can be main- 
tained at the scales of interest since turbulence is subsonic, 
V ^ Cs where v is the turbulent velocity, and so the sound 
crossing time I / Cs (across the correlation length of the tur- 
bulent magnetic fields I) is comparable to or shorter than 
the presumed correlation time of interstellar turbulence 
I /v. Numerical simulations of the multiphase ISM driven 
by supernova explosions indicate statistical pressure equi- 
librium over a wide range of physical conditions in the 
ISM, both between different phases and within a region 
occupied by a single phase (Rosen et al. J,995, Korpi et 
al.HnMI Gazol et al. lMiTji . Anyway, any system that de- 
viates from pressure equilibrium on average must either 
expand or collapse; thus, any statistically steady state of 
the ISM must involve statistical pressure equilibrium. 

Treatment of pressure equilibrium involving magnetic 
field requires some caution since magnetic pressure only 
contributes to force balance across the field lines. However, 
splitting magnetic field into an (isotropic) random part 
and a large scale field, as we do here, admits the reason- 
able assumption that magnetic pressure due to the random 
part is isotropic as long as we consider quantities averaged 
over an intermediate scale which is larger than the turbu- 
lent correlation scale {I = 50-100 pc) but smaller than the 
scale of the regular field (^ Ikpc). 

The peculiar nature of the Lorentz force still remains 
important, e.g. magnetic buoyancy would produce gas mo- 
tions even under perfect pressure balance. However, mag- 
netic buoyancy is important at relatively large scales of the 



order of 1 kpc ('Parker [T979|l . and so can be neglected here. 
Indeed, its characteristic time, equal to the Alfven crossing 
time over the density scale height, cannot be shorter than 
the sound crossing time over the turbulent correlation 
length because both the turbulent scale is smaller than 
the density scale height and the Alfven speed is smaller 
than the speed of sound (again assuming that the turbu- 
lence is subsonic). It is therefore not unreasonable that 
studies of magnetic buoyancy often start with a pressure- 
equilibrium state as an initial condition. Hence, we can 
neglect any effects related to anisotropy of the pressure 
produced by the large-scale magnetic field. 

Pressure equilibrium in the ISM naturally leads to an 
anticorrelation between gas density and magnetic fields 
because regions with enhanced magnetic field must have 
lower density (if only the cooling time is longer than the 
sound crossing time, so that a transient pressure excess 
cannot be removed by cooling). If pressure equilibrium is 
only maintained in the statistical sense, the anticorrelation 
will also be only statistical. 

On the other hand, there are systematic deviations 
from pressure balance in the ISM. For example, the total 
pressure is larger than average in self-gravitating clouds, 
expanding supernova remnants, stellar wind regions, and 
in travelling galactic density waves. Such overpressured 
regions can result in a positive correlation between mag- 
netic field and gas density, leading to enhanced Faraday 
rotation. 

We derive an expression for the Faraday rotation mea- 
sure under the assumption of pressure equilibrium in 
Sect. 13.11 and generalize it in Sect. 13.21 Here we employ 
the assumption of isotropic random magnetic field. As dis- 
cussed in Sect. 14.11 this assumption is only approximate, 
but deviations from statistical isotropy can hardly affect 
significantly the effects discussed in this section. 

3.1. Magneto-ionic medium in pressure balance 

In order to evaluate a correction to RM resulting from the 
anticorrelation of h and rie, we first obtain TIq cLS Oil func- 
tion of h from pressure equilibrium equation. We intro- 
duce Cartesian coordinates (x, y, z) with the a;-axis par- 
allel to B, so that B — (5,0,0), and the line-of-sight 
vector s = {sx, Sy,Sz) directed towards the observer, with 
L the distance from the source of polarized emission to the 
observer. Assuming that cosmic ray pressure is equal to 
the magnetic pressure, we start with pressure equilibrium 
equation, 



where Vg^s is the gas pressure comprising thermal and tur- 
bulent components, the magnetic contribution has been 
doubled to allow for cosmic ray pressure, and V is the 
total pressure. Since Pgas is proportional to the total gas 
density, we write 
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where F = Vg^s/ [Xrit) with nt the total gas density and 
X the degree of ionization. This representation is conve- 
nient if the medium is isothermal, which is a good approx- 
imation since RM is mainly produced in a single phase of 
the ISM, the warm ionized medium (Heiles 1976|. A more 
consistent theory should include fluctuations in the total 
pressure and variable degree of ionization, together with a 
more realistic equation of state (e.g., polytropic). We leave 
these generalizations for future work, and only briefly dis- 
cuss the polytropic equation of state in Sect. 13.21 

Consider magnetic field B consisting of a regular B 
and turbulent parts, and similarly for electron density: 

B = B + b, n^ = N + n, (b) = (n) = ; 

we shall assume that the random field b is isotropic. 
Averaging Eq. Q yields 



—2 



B + (b^) = 47r(-p - NF) 



(5) 




0.0 



2 3 

( 6') /b"" 



Fig. 1. The ratio of RM (solid) and ctrm (dashed) to their 
standard estimates RM*-"-* and crRM(o) given by Eqs H13|l 
and H15|l . respectively. In Eq. H15|l . we have taken / = 0.1 



and B 



Bu 



where we have assumed that V and F do not fluctuate. 
Subtract © from to obtain the following balance equa- 
tion for the fluctuations: 



^2 _ (^2^ +2B b = -AnnF , 



or 



AttF 



2ttF 



-B b 



(6) 



(7) 



This equation shows that electron density is smaller {n < 
0) where magnetic field is stronger, because either is 
larger than average or B and b are similarly directed, 
Bb>0. 

Now we calculate RM identifying the integral in Eq. 
with ensemble average: 



RM/KiL = {B ■ sric) = [B ■ s)N + {{b ■ s)n) 



(8) 



where s is the line-of-sight vector. In order to calculate 
the last term on the right-hand side, we use Eq. to 
obtain 



{nb ■ s) 



1 



2ttF 



{{b-s){B-b)) , 



(9) 



since {{b ■ s) {9 - {9))) = because (6 • s) = and {{b ■ 
s)b'^) = (6'^ cos a) — (with a the random angle between 
b and s) as b is an isotropic random vector. 

The average in Eq. (jS)) can be calculated expanding 
the dot products into Cartesian components: 



{{b-s){B-b)) = ((6,5,-1 
= Bs.Abl) 



bySy 



b,s,)iBb,)) 
Bsy{b^by) + Bs:^{bxbz 



We assume that 6,, by and bz are uncorrelated when eval- 
uated at the same position, so (bxby) = (bxbz) = 0. This 
condition is not restirctive as it is satisfied by any random 
magnetic field b with symmetric probability distributions 
of bx, by and b^ (and for an isotropic b in particular). 
Then 



{{b-s){B-b)) 



ii?ll(&^) 



(10) 



where By = B ■ s = Bsx is the line-of-sight regular mag- 
netic field and |(6^) — {bl) due to the isotropy of b. 

We note that the assumption that the components of 
b are statistically point-wise independent is compatible 
with the solenoidality of magnetic field because the latter 
only requires that magnetic field components are corre- 
lated when taken at distinct positions. In particular, the 
components of a random, isotropic and homogeneous mag- 
netic field must be statistically independent when taken 
at the same position. 

Using Eqs (O and (0, we obtain 



RM = Ki 



B\\NL- 



1 



RM 



(0) 



QttF 



B\\L{b^ 



(11) 



where RM(°) = KiB\\NL is the Faraday rotation measure 
that would arise if magnetic field and electron density fluc- 
tuations were uncorrelated. So, RM has a negative contri- 
bution from the turbulent magnetic field. We can rewrite 
this formula in terms of magnetic field alone assuming that 
the average magnetic pressure is equal to the gas pressure 
(equipartition between magnetic and gas energies). 



+ (&2) ^ 47r(7', 
which yields 



gas/ 



RM RM^"^ 1 - I ^ 




(12) 



(13) 



The equipartition expressed by Eq. I|12|) is assumed to 
hold at scales ^ 1 kpc, so that it docs not contradict the 
assumption of pressure equilibrium (Eq. |0J) that holds at 
small scales, ^ 1 kpc. 



The dependence of RM/RM^°^ on (fe^)/B' is shown in 



Fig. n In the hmiting case {b'^)/B^ ^ 1, we have 
RM = iRM(°) , 
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i.e., B would be underestimated by a factor of 3 if cal- 
culated as B = RM/KiNL as usually done. As dis- 
cussed in Sect. |21 the relative strength of the random 
magnetic field in the Milky Way is {b^y/^/B = 1.3 ± 0.6; 
we adopt {b'^)/B = 2 as a representative value. Then 
Eq. H13() shows that the standard estimate of B is about 
two times too small. In other words, the regular magnetic 
field strength near the Sun consistent with the Faraday ro- 
tation measures of pulsars and extragalactic radio sources 
is i? ~ 3-4 /iG (twice the value usually inferred from RM 
data), in agreement with the estimates from equipartition 
between cosmic rays and magnetic fields. 

Apart from the reduction in the observed RM, quanti- 
fied by Eq. (|13|) . if there were an anticorrelation between 
magnetic field and electron density fluctuations, then the 
observed fluctuations in Faraday rotation measure would 
be reduced as well. When fluctuations in magnetic fleld 
and electron density are uncorrelated, the standard devi- 
ation of RM is given by 

aZ^K^{2(n')(b')Ldf\ (14) 

where d is the common correlation length of magnetic and 
electron fluctuations (see, e.g.. Appendix A in Sokoloff 
et al. 1199 81 . As shown in Appendix ^ under conditions 
where Eq. (|13|l applies, the standard deviation of RM is 
given by 

1/2 



CRM 



40 
9ir 



/ 



,1-./ 



B^ + 



_9_ 

20 ^^2^ 



3 ^11 



1/2 



40 



(62) 



(15) 



where / is the filling factor of thermal electrons, de- 
fined as / = N'^/{n1). (Of course numerical coeffi- 
cients in this formula are model dependent.) Since / ^ 
1 (Berkhuijsen 119991 and references therein), the anti- 
correlation suppresses fluctuations in Faraday rotation, 
crm/cp^i ^ 1- This can be clearly seen in Fig. ^ For 
example, ctrm/ctr^m — 0-4 for {b'^) / B^ ~ 2 and / — 0.1 
(Berkhuijsen 1999). In the limiting case of strong magnetic 

fluctuations, {b"^) / B ^ 1, the right-hand side of Eq. ((T^ 
is approximately equal to f^^^. An important result of 
the reduction of the standard deviation of the Faraday 
rotation measure is that this reduces the depolarization 
effect of internal Faraday dispersion in comparison with 
standard estimates. 

Another implication is that the value of ctrm depends 
on the direction, it is minimum along the regular magnetic 
field and maximum in the directions across it. We stress 
that this anisotropy in ctrm arises even in an isotropic 
random mag netic field. Brown & Taylor ipnTTTjl find that 
the scatter in RM is larger along the Orion arm than across 
it, contrary to our result. As we discuss in Sect. 14.11 it 
is natural to expect that the turbulent magnetic field in 
the ISM is anisotropic in a manner consistent with the 
observations of Brown & Taylor, and that the anisotropy 
in b masks the more subtle effect discussed in this section. 



3.2. Overpressured regions and the polytropic equation 
of state 

Systematic and random deviations from pressure equilib- 
rium are widespread in the ISM, and this affects the re- 
lation between local magnetic field and electron density. 
For example, compression leads to enhancement in both 
magnetic field and gas density, and these variables can be 
correlated. (We note that the anticorrelation discussed in 
Sect. 13.11 is difficult to detect because observational esti- 
mates are biased towards dense regions.) In this section 
we discuss the effect of the positive correlation between b 
and Tie at scales ^ 100 pc, i.e., where the correlation (if 
any) can have a random character; such regions are not 
likely to be widespread in the diffuse ISM. As might be 
expected, this leads to enhanced Faraday rotation. Next 
we consider the polytropic equation of state, where the 
assumption of isothermal medium of Sect. IXTl is relaxed. 

Suppose that the total magnetic field scales with the 
total electron density as 



B^ 
47r 



(16) 



where a and k are constants. The last equality is only valid 
for weak fluctuations, n/N ^ 1. Averaging Eq. (|16|l yields 



B + (6^) = inaN'' , 
b^ - (5^) + 2B b ^ AnaKuN^-^ , 



(17) 
(18) 



which are useful to compare with Eqs (jSJ and © , respec- 
tively. Calculations similar to those in Sect. 13.11 (with F 
replaced by anN'^^^) then yield 



RM = RM(°^ 1 



B 



(19) 



which can be compared with Eq. 113|) . Thus, a positive 
(negative) correlation between B and i.e., k > (k < 
0) results in a positive (negative) bias in the observed 
RM. For example, we get RM = fRM^^^ for k = 1 and 

{b'^)/B 3> 1, i.e., the standard estimate of i? is by a factor 
5/3 too large. 

The corresponding effect on the standard deviation of 
RM is an additional factor k on the right-hand side in 
Eq. fT5|l . 

Now consider a polytropic equation of state, T^gas — 
11,2 Fi, where the factor Fi is similar to F in Sect. 13. ll In 
order to simplify the calculations, we assume that density 
fluctuations are weak, n N,so that n2 — N'' +jnN'*~^ . 
Averaging Eq. Q now yields 



B^ + (6^) ~ 47r(P - 7?i7V'^-iFi) , 

and instead of Eq. 113(1 we obtain Eq. ((19|l with k replaced 
by -7. 
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4. Magnetic field estimates from synchrotron 
intensity 

In this section we discuss how anisotropy of the interstellar 
random magnetic field and the dependence of the cosmic 
ray number density on magnetic field affect the estimates 
of the regular magnetic field from synchrotron intensity. 

4.1. Anisotropic magnetic fields 

Equation (01 applies only to isotropic random magnetic 
fields, but their anisotropy results in stronger polarization 
for a given b/B (LaingJJ81, 2002; Sokoloff et al. 199,8). 
Galactic differential rotation can make the interstellar tur- 
bulent magnetic field anisotropic by extending turbulent 
cells along azimuth so that their azimuthal and radial sizes 
are related via l,p/lr — 1 + AilRAt ~ 1.3 at i? = Rq, 
where Af2 is the angular velocity increment across the 
turbulent cell, R is the galactocentric radius with Rq its 
Solar value, and At ~ 10^ yr is the lifetime of a turbulent 
eddy. This produces anisotropy in the turbulent magnetic 
field, with bcfi/br ~ l,f,/lr- 

A more fundamental reason of anisotropy in b is the 
anisotropic nature of magnetohydrodynamic turbulence 
(see Goldreich & Sridhar 119971 and references therein), 
where turbulent 'cells' are elongated along the mean mag- 
netic field (Goldreich & Sridhar [1995(1 . i.e. roughly in the 
azimuthal direction in the case of Galactic magnetic fields. 

Anisotropy of both types has the largest component 
of the random magnetic field roughly elongated with the 
large-scale Galactic magnetic field b^ 6,- ~ b^. This is 
similar to the anisotropy found by Brown & Taylor H2()()l|l 
in the fluctuations in the Galactic RM. 

The degree of polarization produced by the above 
anisotropic random magnetic field when observed in di- 
rections close to the Galactic centre (where b± consists 
mainly of 6^ and bz with bz — for) is given by [see Eq. (19) 
in Sokoloff et al. [TI 



1) 



Sl + i(62)(a2 + l) 



(20) 



where a = b^/br is a measure of anisotropy. For a = 1, this 
equation reduces to Eq. Q . This estimate also applies to 
directions perpendicular to the spiral arms as long as their 
pitch angle is small. 

The solid line in Fig. [21 shows the fractional polariza- 
tion as given by Eq. 1(20(1 : it is significantly different from 
that given by Eq. ^ for realistic values of parameters. 
The same value of p in fact corresponds to a larger value of 
(6^) /Bj_ when the anisotropy has been taken into account. 
Since the energy density in the random magnetic field is 
limited from above by the kinetic energy of interstellar tur- 
bulence, this implies that values of the regular magnetic 
field strength obtained from Eq. ^ are overestimated. For 
example, the degree of polarization of p/pa = 0.4 is ob- 
tained for {P)/B^ ~ 3.0 according to Eq. and for 




Fig. 2. The fractional polarization with anisotropic ran- 
dom magnetic field as given by Eq. 1(20(1 with a = 1.3 
(solid), under the local equipartition between cosmic rays 
and magnetic field (dashed) according to Eq. 12 1() . and 
from the standard relation 10) (dotted) as functions of 
the ratio {b'^)/B^. The standard relation |j2Jl gives signif- 
icantly overestimated Bj_ for a given fractional polariza- 
tion p and random field strength (6^)^/^ if (b'^)/Bj_ > 2. 

{b^)/B]^ - 2.2 according to Eq. For a fixed (5^), the 
difference in B± is about 1.2. 

Thus, because of the anisotropy of the turbulent mag- 
netic field with b^ ^ br — bz, the regular magnetic field 
obtained from polarized intensity can be significantly over- 
estimated (whereas estimates of the total field remain un- 
affected). In the Milky Way, this effect is strongest in, 
roughly, the azimuthal direction where is dominated by 
&0 and bz and weakest toward the Galactic centre (nearly 
across the spiral arms) where the main contribution to 
b± comes from br and bz- The effect of the anisotropy is 
undoubtedly significant in external galaxies. 

4.2. Local equipartition between cosmic rays and 
magnetic fields 

Arguments similar to those of Sect.|21can be applied to cos- 
mic ray density that may exhibit similar local dependence 
on the magnetic field strength, with similar consequences 
for the magnetic field estimates from the total and polar- 
ized synchrotron intensities. If the energy equipartition or 
pressure balance between cosmic rays and magnetic fields 
were maintained locally (i.e., at any given position), one 
would expect a strong positive correlation between Ucr and 
B (given that Ucr is proportional to the total cosmic ray 
energy density). Sarkar 1(1982(1 argues that such a corre- 
lation, implied by cosmic ray confinement ideas, can be 
enhanced by acceleration of relativistic electrons in com- 
pressed supernova shells. 

In fact, Eq. Q widely used in interpretations of po- 
larized radio emission from the Milky Way and exter- 
nal galaxies is inconsistent with energy equipartition or 
pressure balance between cosmic rays and magnetic fields 
because it is based on the assumption that the number 
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density of cosmic rays is independent of magnetic fields. 
This also applies to Eq. H20I) . If the local equipartition 
is maintained between cosmic rays and magnetic fields, 
so that ricr oc at any position, the total intensity of 
synchrotron emission will strongly depend on the mag- 
netic field since synchrotron emissivity is proportional to 
Uci-B^"' cx B^B^°', where the synchrotron spectral index 
a is close to unity. 

As a result, the equipartition estimate will be biased 
towards regions with stronger field and the degree of po- 
larization will be larger than predicted by Eq. ^ . SokolofF 
et al. [19981 their Eq. (28)] have shown that for n„ oc B^ 
Eq. ^ is replaced by 



lib') 



(21) 

3(&2) + f (&2)2 

where it is assumed for simplicity that the regular mag- 
netic field lies in the sky plane, B = B±, and the ran- 
dom magnetic field is isotropic. The resulting degree of 
polarization is shown in Fig. [3 As expected, Eq. (|21|l 
yields significantly stronger fractional polarization for a 
given ratio (&^)/-Bj_. For example, for p/pQ = 0.4, Eq. (01 

2 _____ 

predicts that {b'^)/B^_ ~ 2.2, whereas Eq. |(2U yields 

— 4.7. For a fixed strength of the random mag- 
netic field, this leads to an overestimate by a factor of 1.4 
in the regular field strength. This difference is very close 
to that between the estimates obtained from synchrotron 
intensities and RM. The equipartition value of the total 
field is also an overestimate. 

A tight, point-wise correlation between ricr and B can 
be an oversimplification. Cosmic ray diffusion can smooth 
variations in ricr at fairly large scales. Then Eq. H21() pro- 
vides an upper estimate of p. Otherwise, one can still use 
Eq. (|21|l . but with (fe^) understood as an average over the 
diffusion length of cosmic-ray electrons. 

There are other uncertainties in the equipartition es- 
timates of field strengths. For example, the standard esti- 
mate relies on the assumption that relativistic protons and 
electrons have the same spectral index in the relevant part 
of the energy spectrum (i.e. around a few GeV for electrons 
responsible for radio synchrotron emission). This is true 
in the Solar neighbourhood (Chapter 9 in Longair 1994). 
However, this is no longer valid if the electrons suffer sig- 
nificant synchrotron or inverse Compton losses, so that 
their energy spectrum is steeper than that of the protons. 
In this case the equipartition or minimum-energy method 
should be applied with care and requires a correction of 
the electron spectrum derived from the synchrotron spec- 
trum. 

5. Discussion 

In order to verify our main results, Eqs. ifT^ and ((T^ . one 
would need to study the relation between RM, P and p 
because Eqs. (|13|l and 1)191) can be rewritten as 



where Eq. H13I) corresponds to k = —1, and Eq. ^ has 
been used with i?^ ~ ^b' and p 'ti Pa- In other words, 
statistical relation between magnetic field and electron 
density fluctuations leads to a correlation of the observ- 
able Faraday rotation measure with the degree of polar- 
ization. In order to rewrite Eq. H22(l in terms of observable 
quantities, we note that -B|| can be related to the observed 
polarized intensity given in Eq. since -B|| and B± are 
related via a geometric factor. Then RM'-"-' can be ex- 
pressed in terms of polarized intensity P, RM*-"-* cx P^/^, 
and Eq. H22() written as 



IRMI 2 

(X 1 — 



3 1 + 3/{2k) po ' 



(23) 



RM ~ RM(") 



^ 3 



Po 



(22) 



where the coefficient of proportionality depends on cosmic 
ray number density, mean thermal electron density, geo- 
metric factors, path length and other variables. Regions 
with pressure balance yield 1 -I- 3/(2k) < 0, and so a 
positive correlation between \RM\/\/P and p, whereas 
overpressured regions should be detectable via an anti- 
correlation between these quantities. Without any inter- 
dependence between magnetic field and electron density, 
|RM|/\/P would be uncorrelated with the degree of po- 
larization p. 

In order to detect the correlation implied by Eq. (|23|) , 
the observational data have to be selected carefully. 
Regions in which synchrotron emission and Faraday rota- 
tion occur together, so that they probe the same magnetic 
field, have to be found. The regular magnetic field must be 
almost uniformly directed throughout the region to reduce 
the variation in the ratio B^\/B±. The positive and nega- 
tive correlations from regions in pressure equilibrium and 
overpressured regions can compensate each other. Thus, 
to detect the correlation produced by pressure balance, 
the region must not be affected by any violent activity 
that could produce overpressure, e.g., expanding super- 
nova shells and Hii regions (Jenkins & Tripp I^UUTIl . The 
region must be very well explored in several wavelength 
ranges in order to have reliable Faraday rotation measures, 
polarized and total synchrotron intensities, and ionized 
and neutral gas densities. Finally, Faraday depolarization 
must not be strong to ensure that P is not affected by 
variations in Faraday depth. 

There are several surveys of the polarized radio back- 
ground in the Milky Way (at A = 21-74 cm - Spoelstra 
[TMl A = 21 cm - Uyamker et al. H^M A = 80-88 cm - 
Haver korn et al. HUM IMI^ and A 21 cm - Gaensler 
et al. 12001(1 . but they cover regions and/or wavelengths 
where Faraday depolarization is strong and so Faraday ro- 
tation and synchrotron emission occur at different depths; 
this makes it difficult to obtain the value of p required. 
Polarization intensity data at smaller wavelengths are 
available (e.g., Duncan et al. 11997(1 . but not the Faraday 
rotation. Further away from the Galactic plane, a few re- 
gions have been observed in radio polarization with high 
resolution, but Faraday rotation data are available only 
in small regions (Reich & Uyamker, in prep.). Further 
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Table 1. Summary of the various effects of random magnetic fields on Faraday rotation measure RM and the degree of 
synchrotron polarization p discussed in this paper. Consequences for the 'standard' estimates of B from the observed 
RM or p are indicated in column 2. Reference to relevant equations is given in the third column; the text should 
be consulted for conditions under which the equations have been obtained; they are briefly summarised in Column 
4, but the effects themselves are more general. Here Uc and Tier are the number densities of thermal and cosmic-ray 
electrons, respectively, with A'^ = (tt-c), and n — Uc — N; B, B and b are the total, regular and random magnetic fields, 
respectively. 



Physical effect 


Consequences 


Equation 
reference 


Assumptions 


TLe-B anticorrelation 


RM reduced, B underestimated, 




pressure balance, isotropic b 








+ magnetic and turbulent energies in equipartition 




(jRM reduced 


m 


as for l|13|l 


n^-B correlation 


RM enhanced, B overestimated 




isotropic b, n/N <^ 1 


Anisotropy of b 


p enhanced, B overestimated 




ricr independent of B 


rici-B correlation 


p enhanced, B and B overestimated 




isotropic b 



radio polarization studies of carefully selected regions in 
our Galaxy are required; Galactic polarization surveys at 
A6 cm are in planning. 

Polarization maps at A < 6 cm are available for several 
external galaxies fBeck 1^0001) . However, the main obstacle 
is the integration along the long line of sight (several kpc) 
and the large beam (several 100 pc) that mixes a range of 
heights above the disc and/or positions within and outside 
the arms. The available resolution of observations is insuf- 
ficient to isolate relatively small regions in the plane of the 
sky where pressure balance can be expected to occur. 

An alternative can be the verification of the effect on 
the standard deviation of RM, described by Eq. H15|l . A 
confirmation of the diminishing effect of pressure equilib- 
rium on fluctuations in the Faraday rotation measure of 
pulsars may have been provided by the results of Mitra 
et al. H2(JU3|I who have shown that the deviations of pul- 
sar Faraday rotation measures from a large-scale model 
distribution are systematically smaller than predicted by 
Eq. (|14|l (we note that the 'estimated arm' in their Fig. 8 
can be misleading as this quantity is obtained by averag- 
ing residuals squared, and so dominated by a few strong 
deviations). The difference appears to be close to a factor 
of two, which is consistent with {b'^)/B ~ 3 (see Fig. QJ. 

The apparent disagreement between the values of the 
regular magnetic field in the Milky Way, obtained from 
Faraday rotation measures and synchrotron intensity, can 
be resolved by allowing for the effects of turbulent mag- 
netic fields, and their correlations with thermal and rel- 
ativistic electron densities. A detailed comparison of the 
various methods to obtain field estimates (radio polariza- 
tion, starlight polarization, extragalactic and pulsar ro- 
tation measures) in a selected field of the Milky Way is 
required to derive reliable estimates of the regular mag- 
netic field in the Milky Way. Regarding external galaxies, 
estimates of the regular magnetic field should be recon- 
sidered with allowance for the anisotropy of the random 



magnetic field and with more consistent implementation 
of the cosmic ray equipartition models. 

An improved estimate of the strength of the regular 
magnetic field in the Milky Way can be important, among 
many other topics, for the studies of the origin of high- 
energy cosmic rays (Kalmykov & Khristiansen I199B|| . 

6. Conclusions 

Our main results are summarized in Table ^ We have 
shown that an anticorrelation between the number den- 
sity of thermal electrons ric and magnetic field strength 
(to which both random and regular magnetic fields con- 
tribute) reduces the Faraday rotation measure RM. We 
have quantified this effect by assuming that the ISM is 
under pressure equilibrium and shown that values of the 
regular magnetic field strength B obtained from relations 
similar to Eq. (j2Jl can be underestimated by a factor of 
about 2. This effect alone can explain the systematic dif- 
ference between the values of B in the Milky Way ob- 
tained from Faraday rotation measures and the degree of 
synchrotron polarization. Another consequence of the an- 
ticorrelation is a suppression of fluctuations in RM. 

On the contrary, a positive correlation between no and 
magnetic field strength results in enhanced RM and over- 
estimated B. This effect can be important along selected 
lines of sight in the Milky Way where the contribution of 
overpressured regions is significant. 

We have also discussed effects that might affect the 
value of the regular magnetic field inferred from the degree 
of polarization of synchrotron emission p. The anisotropy 
of the random magnetic field in the ISM can result from 
the stretching of turbulent cells by Galactic differential 
rotation, in addition to a similar anisotropy inherent to 
magnetohydrodynamic turbulence. The anisotropy of b 
enhances p and can thereby result in B overestimated 
by a few tens of percent. A correlation between the en- 
ergy density of cosmic ray electrons and magnetic field 
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strength also enhances the degree of polarization. We have 
shown that equipartition between cosmic rays and mag- 
netic fields can lead to an overestimate of B by roughly 
the same amount. This implies that the ratio of turbu- 
lent to regular magnetic field strengths obtained from the 
observed degrees of synchrotron polarization should be re- 
vised towards larger values; {b'^)/B ~ 3-4 seems to be a 
plausible estimate. 

There are further effects on the equipartition estimates 
of magnetic field strengths; we briefly discussed the role 
of energy losses of cosmic-ray electrons in Sect. 14.21 

Altogether, the effects discussed in this paper, espe- 
cially the bias in the estimates of B from RM, can reconcile 
the estimates of the regular magnetic field from Faraday 
rotation measures and polarized synchrotron emission. 
However, this requires a more careful interpretation of ob- 
servational data, especially recent observations of diffuse 
synchrotron emission in the Milky Way and new deter- 
minations of Faraday rotation measures of extragalactic 
radio sources and pulsars. 
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Appendix A: The standard deviation of RM 

Here we derive Eq. H15() for the standard deviation of the 
Faraday rotation measure, using the relation = (AT^) — 
(A")^ for the standard deviation of a random variable X. 
Thus, apart from Eq. 0, we have to calculate {{nb • s)^). 
From Eq. Q, we have 

+ {b" - {b^))B ■ b + (B ■ bf] {b-sf . 

The mean values of those terms that contain odd powers of 
b vanish because (i) the components of b are statistically 
point-wise independent (so that, say, {bl-b"^) ~ {bx){bl)) 
and (ii) each of the components has a symmetric proba- 
bility distribution (so, e.g., = 0). Therefore, 

{{b'~{b'))(B.b){b.sr) = o. 

The remaining averages are calculated similarly to the fol- 
lowing: 



{b^ib-sf) 
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where we have taken advantage of the isotropy of b and 
the fact that s is a unit vector. Now we represent in 
terms of its Cartesian components to obtain 

where we used the isotropy of 6, i.e., (fe^) = The 
higher moments of bi are calculated assuming that bi is a 
Gaussian random variable, e.g., (6f) = |(^^)'^ and {bj) = 

Then Eq. H15|l follows after we note that B\\ = Bsx, 
= b\+bI, and 

\inb.sf)-lBl{b^)_ , 

where we have used Eqs (I1U|I and (|12|) . which yields 
F2 = {Vge,s^)/N^ = + (&2))/(47r7V2). We have also 
introduced the volume filling factor of thermal electrons, 
f = Ny{nl)=Ny{N^ + {n^)). 
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